Nonlinear Effects in Pulse Propagation through 
Doppler-Broadened Closed-Loop Atomic Media 
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Nonlinear effects in pulse propagation through a medium consisting of four-level double- A-type 
systems are studied theoretically. We apply three continous-wave driving fields and a pulsed probe 
field such that they form a closed interaction loop. Due to the closed loop and the finite frequency 
width of the probe pulses the multiphoton resonance condition cannot be fulfilled, such that a 
time-dependent analysis is required. By identifying the different underlying physical processes we 
determine the parts of the solution relevant to calculate the linear and nonlinear response of the 
system. We find that the system can exhibit a strong intensity dependent refractive index with small 
absorption over a range of several natural linewidths. For a realistic example we include Doppler 
and pressure broadening and calculate the nonlinear selfphase modulation in a gas cell with Sodium 
vapor and Argon buffer gas. We find that a selfphase modulation of 7r is achieved after a propagation 
of few centimeters through the medium while the absorption in the corresponding spectral range is 
small. 

PACS numbers: 42.50.Gy, 42.65.Sf, 42.65.An, 32.80.Wr 



I. INTRODUCTION 

A main interest in laser driven atomic media is the 
study of their coherence properties. Coherence effects 
like electromagnetically induced transparency (EIT) 
coherent population trapping [!|, lasing without inver- 
sion [H , and others 0, [E] are examples where the opti- 
cal properties of an atomic medium are influenced with 
coherent fields. The interference of different excitation 
channels is the main underlying principle here. A partic- 
ular class of systems in which quantum mechanical inter- 
ference plays a major role are the so-called closed-loop 

systems @, 0, 1 H EE EL HI HI M EE EE El • In 

these systems the laser-driven transitions form a closed 
interaction loop such that photon emission and absorp- 
tion can take place in a cycle. This leads to interference 
of indistinguishable transition pathways between differ- 
ent states. One consequence of this is that it can ren- 
der the system dependent on the relative phase of the 
driving fields. At the same time, however, the inves- 
tigation of closed-loop systems is made difficult by the 
fact that the interfering pathways typically prevent the 
system from reaching a time-independent steady state. 
Such a stationary state in general is only reached when 
the so-called multiphoton resonance condition on the de- 
tunings of the different driving field is fulfilled, which was 
therefore assumed in most previous studies. For general 
laser field detunings, a time-dependent analysis is manda- 
toryiHU. 

Laser driven atomic media are also known to exhibit 
significant nonlinear optical properties IIOL [ill. flp. UA . 

EE El, E3, El, El, HE M HI M, m HE m Mm Ha M. 
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A particular example is the occurrence of an intensity de- 
pendent refractive index, with applications such as beam 
focussing, pulse compression, selfphase- or cross-phase 
modulation or optical switching [IE HE HE HE H3, HI, 
f29l . I30I ] . Here, the connection to coherence properties is 
the following. While an atomic resonance can greatly 
enhance nonlinear effects in atomic media, the accompa- 
nying linear absorption of the same resonance typically 
renders the medium opaque to the probe field. This can 
be overcome by tailoring the response via coherence and 
interference effects. An advantageous situation arises, 
e.g., if the linear absorption vanishes due to destructive 
interference while the nonlinear effect is enhanced by con- 
structive interference. 

Motivated by this, we investigate nonlinear effects in 
pulse propagation through a closed-loop atomic medium. 
In particular, we study a four-level atomic system where 
the four dipole-allowed transitions form a double-A type 
scheme (see Fig. [TJ . Three of the fields are assumed to 
be continous-wave coupling laser fields, while the fourth 
field is a pulsed probe field. We use a time-dependent 
analysis, as the multiphoton resonance condition cannot 
be applied due to the finite frequency spectrum of the 
probe pulses. The medium is modelled as a dilute gas 
vapor including Doppler and pressure broadening and 
an additional buffer gas using realistic parameters. Our 
main observable is the nonlinear index of refraction of 
the medium. We find that our system exhibits a high 
nonlinear index of refraction with small linear and non- 
linear absorption over a spectral range of several natural 
linewidths. In this spectral region of interest, the real 
part of linear and non-linear susceptibility show linear 
dispersion, such that pulse shape distortions are mini- 
mized. For Sodium atoms with Argon buffer gas, we ob- 
tain a nonlinear selfphase modulation of it after 2.9 cm 
of passage through the medium. 

The paper is organized as follows. In the following 
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FIG. 1: (Color online) The four-level atomic system with the 
four dip ole- allowed transitions forming a closed-loop double-A 
type scheme. Three transitions are driven by continuous- wave 
control fields indicated by the solid blue double arrows. The 
fourth transition couples to the pulsed probe field indicated 
by the dashed red double arrow. The coupling strengths are 
given by the Rabi frequencies fljk- The spontaneous decays 
with rates yjk are denoted by the wiggly green lines (j G 
{3,4}, fee {1,2}). 

Sec. Ill Al we present our model. In Sec. Ill Bl we solve 
for the time-dependent long-time limit arising from the 
closed interaction loop in the form of a series. The inter- 
pretation of the series coefficients with respect to their 
physical meaning (Sec. Ill C| will enable us to identify 
the quantities necessary to calculate the linear and non- 
linear susceptibility for the probe field of our system 
(Sec. |HII>|) . Doppler and pressure broadening are dis- 
cussed in Sees. Ill El and III Ft Our results are presented 
in Sec. IIIII both with and without broadening. Finally, 
Sec. IIVI discusses and summarizes our results. 

II. THEORETICAL ANALYSIS 

A. The model 

In this section we present the Hamiltonian for the four- 
level system and the interaction with the coupling fields 
in a suitable interaction picture. We write the field cou- 
pling to transition \j) <-> \k) (j S {3,4}, k e {1,2}) as 

E jk = (e, ]k e-^ + c.c.) , (1) 

with amplitude Ejk, unit polarization vector ejk, and fre- 
quency uijk- For better readability we suppress the space- 
dependence of the fields. The Hamiltonian in dipole and 
rotating-wave approximation reads [4], |5| 

4 

H=^2?iw j A jj 

4 2 , „ 

2 { e ~ tiuJlkt ~ 4 ' lk)A 3k + H.c.} . (2) 

j=3 fe=l 



The energy of level |j) is denoted by hojj and we have 
introduced Rabi frequencies fljk — E^e.^ ■ djk\/h 
with djk being the dipole matrix element of transition 
\j) <-> \k) {j e {3,4}, k £ {1,2}). The complex phase of 
the Rabi frequencies was included into the exponential 
function where tfijk — arg(ejfc • djk). The atomic transi- 
tion operator is defined as Ajk = \j)(k\. 

The canonical approach with a Hamiltonian of the sort 
we have just introduced would be to transform it into 
an interaction picture where the time dependence fully 
vanishes. Unfortunately, this is not possible in our case. 
Due to the closed interaction loop, in general a residual 
time dependence in the Hamiltonian remains. Physically, 
this means that we cannot expect the system to reach a 
true stationary state in the long time limit. The best we 
can do is to use a unitary transformation that gathers 
all the time dependence in a single exponential factor in 
front of the probe field Rabi frequency. In this interaction 
picture we obtain 

Hi =H(A 32 - A 31 )A 22 - hA 31 A 33 
+ h(A 32 - A 31 - A 42 )A A4 

- ^ {Q 31 A 31 + ^32^32 + ^42^42 

+n 41 A 41 e- l(At -^ + H.C.) , (3) 

where the detunings are defined as Ajk = Wjk — {wj — 
We have also defined the so-called multiphoton detuning 
and an equivalent combination of the dipole phases 

A =A 4 i + A 32 - A 31 - A 42 , (4a) 

(j) =041 + ^32 - 031 - 042 • (4b) 

The multiphoton detuning is a typical quantity charac- 
terizing a system with a closed interaction loop. Its sig- 
nificance will become more apparent in Sec. Ill CI 

We now set up the master equation for the atomic den- 
sity matrix g. We include the unitary evolution due to 
the Hamiltonian in the interaction picture and relaxation 
dynamics due to spontaneous decay in Born-Markov ap- 
proximation. The collision induced dynamics will be con- 
sidered in Sec. Ill Ft The unitary evolution is given by the 
Von-Neumann equation and the spontaneous decay can 
be written in Lindblad form j4j. The master equation in 
the interaction picture then reads 

-EEylt^^+H,.}, (5) 

j=3 fe=l 

where g 1 is the density matrix in the interaction picture 
and "/jk is the radiative decay rate of transition \j) <-> |fc). 
For the further analysis we rewrite the master equation 
in a matrix- vector form. Because the trace of the density 
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matrix is conserved we use the corresponding condition 

4 



3=1 



(6) 



to eliminate the diagonal element 5144. Here, 
Qj k = (j'l^lfc). Introducing the vector R = 

(Qii, Bi2i <?13j • ■ ■ i Q43) T containing the remaining fifteen 
elements of the density matrix we find 



dtR + Z = MR, 



(7) 



with an inhomogeneous part E that stems from the elim- 
ination of £144 and a coefficient matrix M. Both E and M 
can be directly derived from the master Eq. ([5]) and con- 
tain the explicit time dependence arising from the time 
dependent Hamiltonian Eq. ([3]). The explicit form of M 
and E is given in the appendix. 



B. Time-Dependent Solution 

To treat the explicit time dependence of the equation 
of motion we first separate E and M into the time in- 
dependent part and the explicitly time dependent part. 
For this, we define 



S =Eq + E_if24ie 



i(At-tf>) 



Eif^ie 



-i(At-tp) 



M =M + M_ifi4ie l(At " 0) + MiCl 41 e 



■i(At-0) 



(8a) 
(8b) 



with time-independent T,j and Mj (j S {0, ±1}). We see 
that under the condition A = the explicit time depen- 
dence vanishes. This is the so-called multiphoton reso- 
nance condition. For fixed coupling field frequencies this 
condition can only be fulfilled for a single probe field de- 
tuning A41 . But we want to investigate probe fields con- 
sisting of pulses with finite temporal length, which due 
to the Fourier relations implies that a whole spectrum 
of probe field frequencies interacts with the medium at 
the same time. Thus, we cannot assume the multiphoton 
resonance condition to be fulfilled Instead, we have 
to solve Eq. ([7]) including the explicit time dependence. 
To do so, we expand R as a power series in ^41, 



R 



OO 

£ 

n=0 



R n £lji 



(9) 



If we assume that the probe field strength is small com- 
pared to the control fields this series will converge. Insert- 
ing Eqs. §5§ and §§§ in Eq. (J7J), we can derive equations of 
motion for the individual coefficients R n . In order ©[OJJ 
we find 

dtR n =MoR n 

+ 6 nA (E_ ie l < At -« + E ie -*( At -«) 

+ (AT-ie^-A + M 1 e^ At ~*A Rn-i ■ (10) 



This is an equation for R n where the coefficient matrix 
Mo is time independent and only the inhomogeneous part 
is time dependent. This time dependence is twofold, first 
again explicitly because of the exponential functions and 
second because of the dependence on R n -i- Thus, we 
make an ansatz for the solution and write R n in a Fourier 
series, 



Rn 



£ 



(11) 



m— — oo 



Projecting on the Fourier basis functions we derive a hi- 
erarchy of time independent equations for the coefficients 
i?l m) . Up to order Oplj] we find 



R Q ° —M 1 E . 



(12a) 



R[ ±1} = (M ± iAty 1 (e±i - M±ii? 0) ) , (12b) 
= - Mq 1 (M^R^ + M x 4" 13 ) , (12c) 



(±2) 
2 

(±1) 



(M ± 2iAl)" 1 M ±1 R\ 
(M ± iAty 1 



(±1) 



(12d) 



x (M ±1 RW + M^R^A , (12e) 

R { ^ 3) = - ( M o ± SiAty 1 M ±1 R { ^ 2) , (12f) 

where 1 is the unit matrix and all other Rn up to this 
order vanish. In general we find that 



oc n 



■im(A.t — <p) 



(13) 



n=0 m • t 
-n+2,. 



Since Fourier coefficients Rn in Eq- 02) only depend 
on Fourier coefficients R^}\ of the next lower order, the 
full solution can be calculated recursively. 



C. Physical Interpretation 

To physically interpret the meaning of the different co- 
efficients we study the influence of the different parts of 
the solution on the probe field. First, we write down the 
expansion series for the relevant probe field coherence in 
the Schrodinger picture Qn using the explicit transfor- 
mation relation connecting the Schrodinger picture with 
our interaction picture. We find 



Q4i =Qli e- 1 ^ 41 *-^ e l(A *-^. 



(14) 



With given as component of the solution for R we 
find 



00 n 



0*1 =£ E K m) 



n=0 m=-r. 

-n + 2,. 
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"41 



x e 



-«[w 4 i+(m-l)A]« i[<t>4i+(m-l)(fl 



(15) 
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where ]i3 refers to the thirteenth component of vec- 
tor Rh\ Thus, coefficient [i?n ]i3 gives a contribution 
at the probe field frequency W41 plus a frequency shift 
of (m — 1)A. The corresponding physical process can 
be identified as follows. A combination of dipole phases 
4> = </>4i — 042+032~ <p3i indicates a full evolution through 
a loop which extends from state |1) to |4) and via |2) and 
1 3) back to state |1). The transition direction is given 
by the sign of the corresponding dipole phase. The evo- 
lution around the interaction loop is also the physical 
reason for the frequency shift A of such a process. Al- 
together, [i2n ]is represents a process with m — 1 loop 
cycles where the sign of m— 1 defines the direction, clock- 
wise for positive or counter-clockwise for negative sign. 
The remaining n — (m — 1) probe transitions can be in- 
terpreted as direct transitions. 



This is not the case for other contributions to ■ F° r 
example, [-Rq ]i3 oscillates at the frequency u> = W41 — A 
and leads to a contribution X ( w = w 3i — w 32 + ^42) 
(four- wave mixing). Here, the resulting frequency is in- 
dependent of LU41. Nevertheless, in principle those pro- 
cesses can influence the result for the linear and third- 
order susceptibility at certain probe field frequencies. For 
example, light can be scattered into the probe field mode 
via different processes. Whether this or similar contri- 
butions change the probe pulse depends on the pulse's 
frequency width compared to the multiphoton detuning 
A and more general also on the propagation direction of 
the probe field relative to the control fields. A definite 
answer to this question requires an analysis of the full 
pulse propagation dynamics through the medium which 
is beyond the scope of this work. 



D. Linear and Non-Linear Susceptibility 



E. Doppler Broadening 



With the above interpretation we can easily identify 
the parts of the solution leading to the linear and nonlin- 
ear susceptibility in the probe field. Because both contri- 
butions should oscillate at the probe field frequency we 
see that m = 1 must be fulfilled in Eq. (|15[) . The order 
of O41 enables one to identify 



oc \r[ 1] 



X (3) (w 4 i) oc 



R 



(i) 
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at O [^Ii] , 
at O [Ojy . 



(16a) 
(16b) 



There is no second order contribution to the suscepti- 
bility as it should be for an isotropic medium [l8l |. By 
comparing the microscopically calculated value for the 
polarization 0, H[ 



P41 =N {d 14 g 41 + c.c.) , 
with the definition of the susceptibility [l£ 



(17) 



=£0- 



E41 



(i) 



£ 4 2 ix (3) 



e 4ie -^i* + c.c. 



we find 



* (1 W) =g^iv 7 4i 

3 „n fa-i , . 3 



R 



13 



(18) 

(19) 
(20) 



with £0 being the permittivity of free space, A41 the wave 
length of the probe field transition, and N the density of 
atoms in the gas. 

We remark that x^H^n) = X^ 3 H W = w 4i~ ^4i+^4i) is 
the lowest order nonlinear contribution at the probe field 
frequency. It leads to an intensity dependent refractive 
index that also depends on oj 4 i and can be different for 
each respective frequency of the probe pulse spectrum. 



A typical experimental setup to investigate the coher- 
ence properties of a laser driven atomic gas would be a gas 
cell with a dilute alkali-atom vapor. For a dilute atomic 
gas theoretical predictions for the linear and nonlinear 
susceptibility can be made on the basis of a single atom 
analysis. This greatly facilitates the theoretical analysis. 
However, in a dilute gas at room temperature or above 
the atoms move at velocities where the frequency shift 
due to Doppler effect cannot be neglected compared to 
the natural line width given by the radiative decay rate 
7. To calculate the Doppler effect for a single field, we as- 
sume a Maxwell-Boltzmann velocity distribution in laser 
pro pag ation direction with a most probable velocity given 

by m 



2k B T 



m 



(21) 



with ks the Boltzmann constant, T the temperature, and 
m the mass of the atom. The non-rclativistic Doppler 
frequency shift is given by 



('-:)• 



(22) 



where uj c g is the shifted frequency seen by the moving 
atom, u> is the lab frame laser frequency, v is the veloc- 
ity of the atom in laser propagation direction, and c is 
the speed of light. The Doppler shift effectively leads to 
an additional detuning Adop with a Gaussian distribu- 
tion m 



/(A Dop ) dA Dop 



1 



dA 



Dop ■ 



(23) 



where k is the wave number. The corresponding line 
width (FWHM) is then given by 



Sw =k\]n(2) 



(24) 
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FIG. 2: (Color online) Real part (solid blue line) and imagi- 
nary part (dashed red line) of the linear susceptibility of the 
probe field. Due to strong control fields Q42 = IOO7 and 
= 5O7 the probe field resonance is split into four dif- 
ferent resonances. Further, Q32 = A31 = A32 = A42 = 0, 
and all spontaneous decay rates 7jfc have been set to 7. The 
susceptibility is plotted in units of 3/8^^4!^. 



To actually calculate the linear and nonlinear suscepti- 
bility for a Doppler broadened medium, for each propa- 
gation direction, we have to add Aoop to the detuning of 
the fields propagating in this direction and then average 
the resulting susceptibility over the velocity distribution 
Eq. 



F. Buffer Gas and Pressure Broadening 

Introducing a buffer gas to the gas cell leads to more 
frequent collisions between the atoms. This has two main 
consequences. First of all it causes pressure broadening. 
For moderate densities, a collision between two atoms 
disturbs the level energies for a short time which results 
in the loss of phase coherence. In a simple approach this 
can be modeled by an additional decay rate j c for the 
coherences. This collisional decay rate consists of a con- 
tribution due to the studied gas itself and a contribution 
due to the buffer gas. Both depend linearly on the re- 
spective density |18j . 

7c =C S N S + C b N b , (25) 

with gas specific constants C s and CV 

A second major effect of a buffer gas is closely con- 
nected to Doppler broadening. Due to the higher density 
the mean free path of a single atom moving in the gas is 
reduced. If it is reduced below the transition wavelength 
an averaging over different velocities during a single emis- 
sion or absorption process can effectively re-narrow a 
Doppler broadened line. This phenomenon is known as 
Dicke narrowing (3lj . 



III. RESULTS 

In principle, Eqs. (fT2j) can be used to calculate analyti- 
cal results for the desired an d X^ 3 ' • But in our situa- 
tion of interest where all four electromagnetic fields, pos- 
sibly all with different detuning, interact with the atom, 
these are usually to lengthy to give any physical insight. 
Therefore, we proceed with a numerical study of the lin- 
ear and nonlinear susceptibility. 

A. Without Doppler Broadening 

Here, our primary goal is to find a set of parameters 
where the intensity dependent refractive index is large 
enough to cause an appreciable amount of nonlinear self- 
phase modulation while the attenuation of a light pulse 
due to absorption is small. To achieve a high non-linear 
index of refraction with low linear and non-linear loss all 
in the same spectral region is challenging because res- 
onances that enhance the nonlinear response typically 
come with strong absorption. Still, we find such a suit- 
able parameter set by manipulating the linear and non- 
linear susceptibility of the probe field as described next. 

We first split the unperturbed resonance of the probe 
field transition by a strong coupling field ^42 and again 
about half as much by the second coupling field £^31 . This 
gives rise to four resonance structures in the linear re- 
sponse, see Fig. [2j 

In this figure, the linear absorption of the resonance 
at A41 ss —257 can be lowered by a small detuning A31, 
which modifies the dressed state populations. Finally, 
optimizing the result with the third coupling we 
can tune one half of the resonance to a small linear and 
nonlinear absorption while still maintaining a substantial 
nonlinear real part. In Fig. [3] it is shown how gradually 
introducing a detuning A31 influences the linear absorp- 
tion, the nonlinear gain, and the real part of the nonlinear 
susceptibility. It decreases the linear absorption and the 
nonlinear gain faster than the real part and thereby im- 
proves their ratio. Interestingly, the imaginary parts of 
the linear and the nonlinear parts of the susceptibility can 
have opposite signs in this spectral region . The linear re- 
sponse induces absorption, while the nonlinear response 
leads to gain. Absorption could in this spectral region 
therefore be reduced even further by a partial cancelling 
of linear absorption and nonlinear gain. However, these 
results are preliminary in the sense, that no effects due 
to Doppler and pressure broadening have been included 
yet. 

B. Including Doppler Broadening 

Using our considerations from Sees. Ill El and IIIFI we 
now want to calculate the linear and nonlinear suscepti- 
bility in a Doppler broadened atomic gas. As a realis- 
tic example we want to assume a Sodium vapor with a 
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FIG. 3: (Color online) Real part (dash-dotted blue line) and imaginary part (solid red line) of the nonlinear susceptibility 
together with the imaginary part of the linear susceptibility (dashed red line). All figures show the resonance around A41 = 
—257. The susceptibility is plotted in units of 3/87r 2 A4 1 A r and for comparability x' 3 ' has been scaled with 3/4£ , | 1 . The 
parameters are A32 = A42 = 0, Q31 = 5O7, Q32 = 347, and Q42 ~ IOO7. The probe field strength is assumed to be one tenth 
of the weakest control field in all cases. The detuning A31 is chosen as (a) A31 = 0, (b) A31 = O.77, (c) A31 = f .57, and (d) 
A31 = I.77. Note the different axis scales in the four subpanels. 



density of N = 1.0 x 10 20 m~ 3 . To reach a vapor pres- 
sure that corresponds to this density the gas cell must 
be heated to a temperature of T = 547.6 K [33| . At this 
temperature the Doppler linewidth is 8uj = 2ttx1.78 GHz 
which is very broad compared to the natural linewidth of 
the Sodium Di transition of 7 = 2irx 9.76 MHz. In a pure 
Sodium vapor the spectral features we found in Sec. lIII Al 
would be averaged out by the Doppler effect. But if we 
introduce a buffer gas strong pressure broadening can 
preserve them. For Argon and Sodium, the gas param- 
eters in Eq. (j25|) are given by C s = 1.50 x 10~ 13 m 3 s _1 
and C b = 2.53 x 10 -15 m 3 s -1 [Hj|. We want to as- 
sume a collision-induced coherence loss rate of 7 C = 
1.0 GHz which corresponds to a buffer gas density of 
Nb = 3.95 x 10 23 m~ 3 . At such a density the mean free 
path is of order A = 10~ 5 m. This is much larger than 
the transition wavelength A = 589.2 x 10 -9 m such that 
the limit of Dicke narrowing is not reached. 

We now try to recover results similar to the unbroad- 
ened case shown in Fig. [31 Because of the strong broad- 
ening we have to apply correspondingly stronger control 
fields. For fl 42 = 60.0 GHz and n 31 = 30.0 GHz, we 
find the resonance studied in the unbroadened case at 
around A41 = —15.0 GHz. The third control field is set 
to f2 32 = 25.0 GHz and the detuning to A31 = 1.6 GHz. 
For the Doppler averaging we have assumed all fields to 
be co-propagating. The different subpanels in Fig. 3] cor- 
respond to different Doppler linewidths, and thus via 



Eq. (HM)) to different temperatures. In Fig. 0Ja), the 
Doppler linewidth is chosen below the natural linewidth 
of the probe transition, and as expected we finds results 
that are similar in shape to the unbroadened case (see 
Fig-Eld)). Differences are mainly due to pressure broad- 
ening. Gradually increasing the Doppler linewidth up to 
the full Doppler width expected for the gas parameters 
discussed above in subfigure (d), we find that while the 
shapes of the different curves change, our main result 
of high nonlinear index of refraction with small linear 
and non-linear absorption persists with Doppler broad- 
ening. Also in the broadened partial cancelling 
of linear absorption and nonlinear gain could be possible. 
Note that since the averaging process affects not only the 
probe field detuning but all four detunings at the same 
time the results cannot be explained in terms of a simple 
smoothing of the curves without Doppler effect. 

We also considered different laser geometries, such as 
control fields propagating perpendicular to the probe 
field, or one or two control field propagating in oppo- 
site directions, and found the co-propagating case to be 
the most advantageous one. This is similar to the case 
of Doppler broadening in typical electromagnetically in- 
duced transparency setups where co-propagating lasers 
typically are preferable. 

We finally use our results at probe field detuning A = 
— 17.8 GHz to calculate the required optical length for 
a nonlinear selfphase modulation of tt. This probe field 
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FIG. 4: (Color online) Real part (dash-dotted blue line) and imaginary part (solid red line) of the nonlinear susceptibility 
together with the real part (blue dotted line) and the imaginary part of the linear susceptibility (dashed red line) at the 
resonance around A41 = —15.0 GHz. The control fields have Rabi frequencies £742 = 60 GHz, = 30 GHz, Q32 = 25 GHz, 
and the detunings are A31 = 1.6 GHz, A32 = A42 = 0. The medium parameters described in the main text correspond 
to Sodium as the active medium with Argon as a buffer gas. The four different plots show Doppler averaged results with 
a Doppler linewidth of (a) below the natural linewidth, (b) 50%, (c) 90%, and (d) 100% of the full Doppler linewidth of 
Slj = 2tt x 1.78 GHz. 



frequency is indicated by the vertical blue dotted line in 
Fig. [|Jd). The nonlinear selfphase modulation is given 
by Gf 

A$ m =n 2 IkL, (26) 

with 712 the intensity dependent refractive index, / the 
probe field intensity, k the wavevector, and L the optical 
length. We assume a probe field strength one tenth of 
the smallest control field and find 

L„ =2.9 cm . (27) 

From Fig. [4] (d) we see that the magnitude of the imag- 
inary parts of the linear and nonlinear susceptibility are 
more than one order of magnitude smaller. Therefore, 
the equivalent characteristic length scale is more than 
one order of magnitude larger. Furthermore, both parts 
give rise to small gain rather than absorption. 

Thus, our results show, that in a certain spectral region 
a nonlinear selfphase modulation of 7r can be achieved on 
a realistic laboratory lengthscale. Since the real part of 
both the linear and the nonlinear susceptibility display 
approximately linear dispersion in the spectral region of 
interest, pulse shape distortions can be expected to be 
small. Interestingly, the real part of the linear suscep- 
tibility has a negative slope in the considered frequency 
region, in contrast to a positive slope typically found in 
an electromagnetically induced transparency window. 



IV. CONCLUSION 



We have studied nonlinear effects in pulse propaga- 
tion through a laser-driven medium where the applied 
fields form a closed interaction loop. Such loop systems 
in general only allow for a time- independent treatment at 
a single probe field frequency, where the so-called multi- 
photon resonance condition is fulfilled. As a probe field 
pulse has a finite frequency width, this condition which 
allows for a straightforward theoretical treatment could 
not be applied. Instead, we treated the time-dependent 
problem by turning it into a hierarchy of equations that 
describe the various physical processes occurring in the 
medium. We have included Doppler and pressure broad- 
ening as well as a buffer gas in our analysis and have used 
realistic parameters for a medium consisting of Sodium 
vapor. We could show that the studied system can ex- 
hibit a high non-linear refractive index with small ab- 
sorption or gain over a spectral range of several natural 
line widths. For the chosen parameters, both the linear 
and the non-linear susceptibilities show near-linear dis- 
persion such that pulse shape distortions are minimized, 
and the slope of the linear dispersion is negative. A non- 
linear selfphase modulation of n is obtained after 2.9 cm 
propagation through the medium. 
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APPENDIX: COEFFICIENTS MATRIX 

The explicit form of the coefficient matrix M and the 
inhomogeneous part E can be derived from Eq. frj). Here, 
we list all nonzero elements Mj ^ and Ej , which are given 
by 

Mi tl = M 1)6 = M 6 ,6 = 2 M ii>n 

= Ei = ^6 

= -7r, 

M^g = M x , 9 = M 2 ,w = M 3 ,4 

= M 4 ,i2 = M* 7 = M 9 ,n = M* 315 

M* 3 = M 5>9 = M* 7 - M 6 ,io 

= M 7) n = M 8 ,i2 = Mr ,n = M* 415 

= 2^32 , 

M* 4 = Mi,i3 - M a ,i4 - M 3) i B 

= -M^ = M 4 * 6 - M 4 * u = M* 8 

= M 9,12 = ^ f 13,l = Mi3, 8 = Mi3, U 

= E4 = Ei 3 

= 1 -i(A*-« 

2 

M 2 * 4 = M 5 ,i3 = M* s = M 6 , 14 



= M 7 ,i8 = M 8 * ! = ~M* 6 = M* n 

= ^10,12 = -^14,1 = 2^14,6 = Afl4 : ll 

= E 8 = E14 

= 2^42 , 

M 3 , 3 = M* 9 

= -7 r - iA 31 , 
M 4 , 4 = M* 313 

= -7 r - i(A 3 i + A 42 - A 32 ) , 
M 7 , 7 = M* aAQ 

= -~f r - iA 32 , 
Af 7i8 = Mio.io 

= -7 r - iA 32 , 
A/12,12 = M* 5)15 

= -27 r - i(A 42 - A 32 ) , 
M 2i2 = M 5 * 5 

= -i(A 3i - A 32 ) , 

M 6i4 = Mn,4 = M 6 ,13 = 13 

= Mi,b = Mn >8 = Mi, 14 = Mu M 
= 0, 



where Mj fc = Mkj holds if not noted otherwise and by 
M* k we indicate the complex conjugate of M^/.. 
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